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CHAPTER  I 


INTRODUCTION 

In  order  to  investigate  the  scattering  effects  associated  with 
antennas  in  the  presence  of  scattering  structures,  two  computer  codes 
have  been  developed:  The  Aircraft  Antenna  Code  [ll  and  The  Basic 
Scattering  Code  [2l.  The  scattering  structures  associated  with  these 
codes  are  modeled  by  finite  flat  plates  and  cylinders,  e.g.,  the  wings 
or  tail  structures  of  an  aircraft  can  be  modeled  by  plates  and  the 
fuselage  by  a  cylinder.  Both  codes  use  geometrical  theory  of  dif¬ 
fraction  (GTD)  to  analyze  the  scattering  effects,  which  are  most  use¬ 
ful  for  low  gain  antennas  although  they  can  simulate  large  antennas 
using  the  superposition  principle.  Thus  scattered  fields  can  hp 
calculated  by  superimposing  the  contribution  from  each  element. 
Unfortunately,  this  approach  results  in  large  computer  run  times  for 
antennas  that  have  to  be  subdivided  into  many  smaller  elements. 

Two  approaches  for  overcoming  this  limitation  are  investigated 
in  this  thesis:  a)  the  large  subaperture  technique  and  b)  ttie  two- 
point  GTU  source  field  method.  Both  approaches  use  the  GTD  methou 
to  calculate  the  scattering  from  structures  in  a  manner  similar  to 
that  used  in  the  present  codes[l,2l;  however,  the  fields  of  the  an¬ 
tennas  or  sources  are  treated  in  a  way  to  permit  more  efficient  cal¬ 
culations  of  these  approaches.  Both  approaches  are  useful  and  com- 


plement  each  other.  For  the  purpose  of  review,  the  basic  GTD  theory 
is  presented  briefly  in  Chapter  II  and  a  two-dimensional  model  of 
the  existing  computer  codes  is  discussed  in  Chapter  III. 


The  first  approach  is  discussed  in  Chapter  IV  which  uses  large 
subapertures  to  represent  the  large  antenna  aperture  distribution. 
Basically,  it  calculates  the  antenna  fields  in  the  same  way  as  the 
aperture  integration  (AI)  method  except  that  elements  are  grouped 
into  large  subapertures.  This  method  works  well  in  most  regions  if 
the  large  subapertures  are  no  larger  than  one  wavelength.  Otherwise, 
it  only  works  well  in  some  regions,  especially  near  the  main  beam. 

In  the  second  approach,  which  is  discussed  in  Chapter  V,  GTD 
provides  a  very  efficient  method  because  the  antenna  fields  can  be 
modeled  by  a  2-point  edge  diffraction  technique  in  the  regions  suf¬ 
ficiently  removed  from  the  main  beam.  However,  GTD  can  not  be  used 
to  calculate  the  antenna  fields  in  the  main  beam  since  one  loses  some 
information  about  the  aperture  distribution  between  the  two  edges 
of  the  antenna. 

In  the  two-point  GTD  source  field  approach,  the  sources  of  the 
incident  fields  on  the  strip  scatterer  are  located  at  each  edge  of 
the  antenna.  When  a  strip  scatterer  is  near  the  projected  aperture 
of  the  antenna,  the  incident  field  from  the  source  at  the  nearest 
edge  has  a  large  slope.  Thus  ordinary  edge  diffractions  from  the 
scatterer  become  insufficient  to  get  a  good  result.  So,  slope  dif¬ 
fraction  from  the  scatterer  was  added  to  partially  solve  the  problem 
as  discussed  in  Chapter  VI.  However,  when  a  strip  is  very  close  to 
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a  shadow  boundary  of  the  projected  aperture,  the  incident  field  from 
the  edge  source  is  too  complicated  to  model  by  a  uniform  wave  plus 
a  slope  wave.  Consequently,  the  large  subaperture  approach  must  be 
used  for  these  cases. 

The  results  and  discussions  are  presented  in  Chapter  VII  to  show 
the  efficiency  and  limitations  of  the  two  approaches.  Finally,  two 
specific  examples  are  demonstrated  in  Chapter  VIII  so  that  the  reader 
can  relate  this  study  to  practical  application. 

In  this  thesis,  only  the  two-dimensional  (2-D)  case  is  investi¬ 
gated  in  that  it  provides  insight  for  the  problem  without  adding  the 
unnecessary  complexities  associated  with  the  three-dimensional  (3- 
D)  case.  For  the  3-0  case,  some  more  complicated  factors  such  as 
geometry,  polarization,  and  corner  diffraction  must  be  considered. 
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CHAPTER  II 


BASIC  GTD  THEORY 

The  basic  GTD  solutions  of  the  two-dimensional  case  are  reviewed 
in  this  chapter,  and  these  solutions  are  applied  to  solve  fields  in 
the  following  chapters.  When  a  wedge  is  illuminated  by  a  line  source, 
GTD  is  used  to  analyze  the  radiation  from  the  source  and  the  scat¬ 
tering  from  the  wedge.  For  the  geometry  shown  in  Figure  1,  the  total 
field  caused  by  the  line  source  in  the  presence  of  the  wedge  is  ex¬ 
pressed  by 

=  Es  +  Er  +  Ed  (2-1) 

The  source  field  Es  is  the  electric  field  radiated  directly  from  the 
source,  the  reflected  field  Er  is  the  electric  field  reflected  from 
the  surface  of  the  wedge,  and  the  superposition  of  these  two  fields 
is  called  geometrical-optics  field.  The  diffracted  field  Ed  from 
the  edge  of  the  wedge  includes  ordinary  edge  diffraction  Ed,  and  slope 

i 

diffraction  E<j.  For  most  cases,  first  order  diffraction  is  enough 
to  solve  the  scattering  problem,  but  in  some  cases,  slope  diffraction 
must  be  added  to  the  solutions. 
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OBSERVATION 

POINT 


Figure  1.  Geometry  for  two-dimensional  wedge  diffraction. 


A.  Geometrical-Optics  Field 

As  shown  in  Figure  1,  the  incident  shadow  boundary  (1SB)  is  de¬ 
termined  by  the  locations  of  the  line  source  and  the  edge,  and  re- 
I  lection  shadow-boundary  (RSB)  is  determined  by  the  image  source  and 
the  edge.  Once  these  boundaries  have  been  defined,  the  source  field 
Es  can  be  expressed  by 


(2-2) 


-jkps 

'Fs . 


for  0  <  4>  <  100  +  <t>* 


elsewhere 


That  is,  the  source  field  exists  in  regions  I  and  II  as  shown  in  Fig¬ 
ure  1.  Similarly,  the  reflected  field  Er  can  be  expressed  by 
-jkPr 

+  ~~ -  for  0  <  <|>  <  180°  -  4)' 


(2-3) 


elsewhere 


which  exists  in  region  I  only.  Note  that  pr  is  the  distance  from 
the  image  source  to  observation  point  and  one  uses  the  minus  sign 
for  the  TE  case  and  plus  sign  for  the  TM  case. 

B.  Ordinary  Edge  Diffraction 


1.  The  wedge 


The  ordinary  diffracted  field  at  the  observation  point  (s,4>) 


is  given  by  [3] 


Eo  =  eI(QE)  -  Ds,h(L;4>><*>’ »n)  7^- 


(2-4) 


where  E^(Q^)  is  the  incident  field  from  the  line  source  at  the  edge 
Q^,  which  can  be  expressed  by 


E1 (Qe)  = 


e-jks' 

yi7 


(2-5) 


Note  that  s'  and  s  are  the  distances  between  the  edge  and  the  source 
and  the  observation  point,  respectively.  The  phase  factor  is  e-^5 
and  l/JT  is  the  2-0  spreading  factor.  The  diffraction  coefficients 
Ds  and  Dh  are  given  by 


D$(L; <j>,<)>'  ,n)  =  Dj (L;<{>-<|>' ,n)  +  D ^ ( L ; 4>+4> '  »n)  (2-6) 
h 

where  the  minus  (-)  sign  applies  for  Ds,  which  corresponds  to  the  TE 
case,  and  plus  (+)  sign  applies  for  D^,  which  corresponds  to  the  TM 
case.  Note  that  Dj  is  given  by 


IT 


+  .  wA 
where  0=  4>-4>  and  n  =  2  -  — , 

with  wedge  angle  WA  given  in  radians.  The  term  a±(0)  is  a  measure 
of  the  angular  separation  between  the  observation  point  and  the  in¬ 
cident  or  reflection  shadow  boundary  and  is  expressed  by 
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(2-8) 


a'(g)  =  2  cos2(^Lj 
+ 

where  N~  are  integers  which  most  nearly  satisfy  the  equations 

2nnN+  -  8  =  n  ,  (2-9) 


and 


2imN~  -  0  =  -  it 


(2-30) 


The  distance  parameter  L  is  given  by 


L  = 


ss* 

s+s1 


(2-11) 


and  the  transition  function  F(x)  is  basically  a  Fresnel  integral, 
which  is  given  by 


F(x)  =  2j Jx  eJK  f 


:2 


dr 


(2-12) 


Jx 


where  X  =  kLa  (8) 


If  the  line  source  is  far  away  from  the  edge,  that  is,  in  the  case 
of  s'-**®,  the  incident  field  becomes  a  plane  wave.  In  that  case,  the 
distance  parameter  L-*-s. 
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2.  The  half  plane 


When  n=2,  the  wedge  becomes  a  half  plane  and  the  diffraction 
coefficients  D$  and  Dh  can  be  simplified  to 

D  (LiM1)  =  DjU;*-*1)  +  .  (2-13) 

In  which  Dj  are  expressed  by 
_j  JI 

0i"-;b)  *  fpjnr  <?-14> 


and  a(g)  is  simply  given  by 

a(B)  =  2  cos2  |  (2-15) 

C.  Slope  Diffraction 

When  there  is  zero  incident  field  on  the  edge  or  the  slope  of 
the  incident  wave  is  large,  slope  diffraction  becomes  significant. 

This  section  derives  the  equations  for  the  slope  diffraction,  which 
has  a  form  similar  to  ordinary  edge  diffraction  except  that  the  slope 
diffraction  coefficients  9DS  3<{> '  and  the  slope  of  the  incident  field 
at  the  edge  9E^(Q^)/3n  are  used. 
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Slope  diffraction  for  a  cylindrical  wave  can  be  derived  from 
edge  diffraction  by  considering  a  line  dipole  source  composed  of  two 
line  sources  [4]  as  shown  in  Figure  2.  The  field  of  this  dipole  source 


is  given  by 


E 


1 


I  • 


(2-16) 


where  s+  and  s  are  the  respective  distances  from  each  line  source 
to  the  field  point.  For  small  spacings  £<< s * ,  one  obtains 


>+=  s'+  i  Hsi 


(2-17) 


where  the  angle  i|>  is  shown  in  Figure  2.  Thus,  the  source  field  can 
be  expressed  as 


(2-18) 


Figure  2. 


A  line  dipole  source  for  slope  diffraction. 
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T 


I  nr  .1  slope  <1  i  1 1  r.ir  I.  ion  source  k(t>0,  and 

-jks 1 

L1  :  jk  IH  s  i nrp  e -  .  (2-19) 

JT7 


The  slope  diffraction  field  can  be  derived  by  superimposing  the 
diffracted  fields  from  the  two  line  sources  as  shown  in  Figure  3. 
Thus  using  Equations  (2-4),  one  finds 


-jks’ 


s' 


+  ”  Ds,h(L’<t>’4> '  " 


•jks 


s 

(2-20) 


Since  A4>  '-*-0  as  l«s',  Equation  (2-20)  can  be  expressed  as 


e- .  ,Ds".  ...  e 

Ls  i—  3<{>  '  W 


-jks 


7s' 


(2-21) 


Let  n  be  the  distance  normal  to  the  edge  and  normal  to  the  in¬ 
cident  ray  from  the  source  as  shown  in  Figure  4.  Note  that  n+s '  •  (tt-!|') 
as  .  Then  the  slope  of  the  incident  field  at  the  edge  can  be 
expressed  in  terms  of  the  distance  n  as  follows: 


IE1!  .  -1  aL1  I  -  jkU  e'Jks' 

3n  \  r  K  s'  fs7 


Since  l  -  s'  A4>',  Equation  (2-22)  can  be  expressed  as 


I  A4>* 


e-jk5>' 


(2-22) 


(2-23) 
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POINT 


'X 
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Figure  3.  2-D  slope  diffraction  for  a  wedge. 


Figure  4.  Normal  vector  n  of  the  incident  field  for 
slope  diffraction. 


Combining  Equations  (2-21)  and  (2-23),  one  gets 


E(l  =  *  -2L.| 

s  ji<  an 

UE 


(2-24) 


The  slope  diffraction  coefficient  aD$  ^/’3<|)'  can  be  expressed  by 


3D 


(2-25) 


where 


DpiU;e) 


30,(L;(!) 
3(j) 1 


(2-26) 


Note  that  the  plus  sign  is  used  for  the  TE  case  and  minus  sign  for 
the  TM  case,  and  Dj  is  defined  in  Equation  (2-14). 
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CHAPTER  III 


APERTURE  INTEGRATION  SOURCE  AND  GTD  SCATTERING 

The  original  method  used  to  solve  a  2-D  scattering  problem  is 
studied  in  this  chapter.  The  total  scattered  fields  are  calculated 
in  the  same  manner  as  the  existing  codes  [  1,2]  for  the  purpose  of 
reference  and  comparison.  Figure  5  shows  the  2-D  model  which  is 
analyzed  in  the  following  chapters.  A  source  antenna  of  width  L  is 
located  between  edge  points  Qsi(xspysi)  and  QS2^Xs2’-ys2^  and  a  str’P 
with  edge  points  Qcl ( xcl »yci )  and  rePresents  a  perfectly 

conducting  scatterer.  Then  the  total  field  of  the  antenna  with  the 
strip  pattern  present  is  composed  of  the  source,  the  reflected,  and 
the  diffracted  fields,  which  are  calculated  separately. 

The  fields  scattered  from  the  strip  are  calculated  for  each  ele¬ 
ment,  that  is,  each  element  of  the  antenna  aperture  distribution  is 
taken  as  a  source  for  the  diffracted  fields  as  shown  in  Figure  6. 

Then  the  contributions  from  each  element  are  superimposed  to  obtain 
the  total  field.  This  method  gives  accurate  results,  but  is  not  ef¬ 
ficient  for  large  antennas.  Two  methods  for  improving  the  efficiency 
are  discussed  in  Chapters  IV  and  V. 
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A.  Source  Field  Without  the  Scatterer 


Consider  the  fields  of  the  antenna  without  the  scatterer  first. 
Using  the  aperture  integration,  the  antenna  is  divided  into  N  small 
elements.  For  the  nth  element  dyp  located  at  Qsn(0>ysn)  as  shown 

c 

in  Figure  6,  the  field  contribution  dEn  is  given  by 


A(y 


sn 


dyn 


„  e 
cos  -g 


(3-i) 


0 

where  A(y)  is  the  aperture  distribution  and  a  cos  element  pattern 
gives  the  most  appropriate  modeling  for  an  edge  as  derived  in  Appendix 
A.  This  approach  provides  excellent  overlapping  of  GTD  and  AI  for 
the  analysis  of  reflector  antennas  [5,6].  Since  s'  =  R  -  ysn  sine 
for  the  far  field,  one  finds 


dEn  =  A(ysn)  ejkysnsine  dyn  .  cos  | 


,-jkR 


(3-2) 


So  the  total  far  field  pattern  function  F(e)  is  obtained  by  super¬ 
posing  all  the  contributions  from  the  N  elements,  such  that 

F(0)  =  f  A(ys_)ejkysnsin9  Ay  cos  |  (3-3) 

n=l  d 

where  Ay  =  L/N 


and  ysn  =  *  ^  +  2*  +  (n-l)Ay. 
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B.  Total  Field  with  the  Strip  Scatterer  Present 

Next,  consider  the  antenna  in  the  presence  of  the  strip  scat- 

c 

terer  as  shown  in  Figure  7.  Let  0nj  represent  the  angle  and  snj  rep¬ 
resent  the  distance  from  Q$r)  to  Qcj,  respectively,  where  n  =  1,  2, 

...,  N  and  j=l,2.  and  0^  correspond  to  the  shadow  boundaries 

of  the  nth  element  at  Q  .  Then  the  source,  the  reflected  and  the 
diffracted  fields  can  be  analyzed  as  follows. 


Figure  7.  The  source  field  and  its  shadow  boundaries 
of  the  nth  element. 

1 .  Source  field  with  the  scatterer  present 
Since  some  of  the  source  fields  are  blocked  by  the  strip  scat¬ 
terer,  it  can  be  seen  from  Figure  7  that  the  pattern  function  of  the 
source  field  from  element  Qsn  can  be  represented  by 
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Fn(  e)  -  < 


for  e*!  <  e  <  6^2  »  and 


jkysnsin0 
^A(y$n)e  Ay  elsewhere 


(3-4) 


fly  suporpos  it  ion,  the  total  source  field  pattern  is  given  by 


Fs(e) 


N  c 

l.  ■>) 

n=l 


(3-5) 


Note  that  each  source  element  has  its  corresponding  shadow  boundary 
angles  at  and  e^* 

2.  Reflected  field  from  the  strip 

At  the  same  time,  there  are  some  fields  reflected  by  the  scat¬ 
tered  The  image  method  is  used  to  calculate  the  reflected  fields. 

As  shown  in  Figure  8,  the  nth  image  element  is  located  at  Qin(x{n>yin) 
and  its  contribution  to  the  reflected  field  is  given  by 

J'P  i  n  R  R 

CR,.  +A<*sn>  e  ,orV<  e  <  e„i  .  and 

V6  '  M 


(3-6) 

(0  elsewhere 


whore  the  sign  is  defined  in  Equation  (2-3)  and  jn  =  k(xIn  cose+yIn 

sinel  refers  the  phase  to  the  origin  (0,0).  The  shadow  boundary  angles 
R  R 

are  0(1j  and  0^  of  the  image  point  QIr|.  Then, the  total  pattern  of 
the  reflected  field  is  given  by 


18 


fr(0)  =  I  F*(e) 

n-1  n 

3.  Diffracted  field  from  the  edges  of  the  scatterer 

Next,  the  diffracted  fields  from  the  two  edges  of  the  scatterer 
are  evaluated.  For  the  nth  source  element  Qsn,  there  are  two  dif¬ 
fracted  fields,  E^(0)  and  E^e),  from  the  two  edges  of  the  scat¬ 
terer,  respectively,  as  shown  in  Figure  9.  The  far  field  pattern 
F^j(e)  of  the  diffracted  field  from  Qcj.  as  caused  by  the  nth  source 
element  is  given  by 


Figure  9.  The  geometry  for  the  diffracted  fields. 


¥ 


I 


where  the  incident  field  E1  .  from  0  to  0  .  is  qiven  by 

nj  sn  ^cj  3  J 


_jksni 

I  /  r.  \  -  A  /  \  e  0 

— —  cos  ^ 


£nj<0cj)  *  "(ysn) 


(3-9) 


I  s„ . 

nj 


and  D  •  is  the  diffraction  coefficient  which  is  given  by 

*  *  J 


D  .  =  Dr(L  .Hb'.)  +  ML  -  ;4>  .+4>' •) 
nj  P  nj,vnj^nj;  I1  nj’^nj  Ynj' 


(3-10) 


where 


D j ( L ; 4>  — 4> 1 )  was  defined  in  Equation  (2-14). 


Note  that  one  should  use  the  sign  for  the  TE  case,  and  the  "+" 
sign  for  the  TM  case.  Further, 


-nj 


Snj*  Sj 
Snj+Sj 


s  •  as  s  .-hx> 
bnj 


and 


sj  =R-(xcjCos0+ycjsin0) 

So  the  total  pattern  function  associated  with  all  the  first  order 
edge  diffracted  fields  is  given  by  summing  all  the  diffracted  fields 
from  the  N  elements,  such  that 
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(3-11) 


Fd(6) 


?  f 

n=l  j=l 


Enj^cj^ 


dj 


whore 


vi)dj  =  k(xcjCOS(3+ycjSin0) . 


(3-12) 


4.  Total  far  field 

Recall  that  the  far  field  pattern  of  the  antenna  with  the  strip 
scatterer  present  is  given  by  superimposing  the  source,  reflected, 
and  diffracted  fields  as  stated  in  equation  (2-1).  Thus, 

Ftot(e)  =  Fs(0)  +  FR(e)  +  Fd(e). 
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CHAPTER  IV 


LARGE  SUBAPERTURE  TECHNIQUE 


When  an  antenna  or  source  aperture  is  quite  large,  it  takes  a 
lot  of  computer  run  time  if  one  uses  the  method  discussed  in  the  last 
chapter.  Most  of  the  time  is  spent  evaluating  the  diffracted  fields 
for  each  little  source  element.  For  better  efficiency,  the  large  sub¬ 
aperture  technique  groups  the  source  into  equal  size  subapertures. 

The  contribution  of  each  subaperture  to  the  diffracted  field  is  ob¬ 
tained  by  using  the  subaperture  as  a  single  source.  In  other  words, 
its  incident  field  on  the  edges  of  the  scatterer  is  the  sum  of  the 
contributions  from  all  the  little  elements  within  the  larger  limits 
of  the  subaperture;  these  contributions  are  calculated  in  the  same 
way  as  done  in  the  previous  chapter.  Then,  the  mid-point  of  each 
subaperture  is  used  to  evaluate  the  diffraction  coefficient  for  that 
subaperture.  This  approximation  works  well  as  long  as  the  subapertures 
are  not  too  large. 

The  source  and  the  reflected  fields  are  grouped  with  the  dif¬ 
fracted  fields  so  that  the  total  field  is  continuous  across  the  shadow 
boundaries.  This  will  be  shown  by  computer  plots  given  in  Chapter 
VII.  The  analyses  of  the  source,  reflected,  and  diffracted  fields 
are  similar  to  those  in  Chapter  III  and  given  in  the  following  sec¬ 
tions. 
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A.  Source  Field  with  the  Scatterer  Present 


In  this  approach,  each  subaperture  which  includes  several  ele¬ 
ments  is  considered  as  an  individual  source.  Let  M  be  the  number 
of  little  radiators  within  the  subaperture  limits.  The  source  field 
pattern  of  the  kth  subaperture  is,  then,  given  by 

[  0  for  0^  <  0  <  0k2»  and 


(4-1) 


M  jkykfnsin0 

l  Mykm)e  *Ay  elsewhere, 

m=l 

V 


where  ykm  represents  the  location  of  the  mth  element  in  the  kth  sub¬ 
aperture,  and  0kj,  0^  are  the  shadow  boundary  angles  which  are  de¬ 
termined  by  the  mid-point  (0,yk)  of  the  kth  subaperture  as  shown  in 
Figure  10. 

/ 

/ 

/ 

/ 


Figure  10.  The  source  field  and  its  shadow  boundaries 
for  the  kth  subaperture. 


Then,  the  total  source  field  is  the  superposition  of  the  source 


fields  of  all  subapertures  as  given  by 

FS{0)  =  l  F?(6)  , 
k=l  K 


(4-2) 


where  K  is  the  number  of  subapertures. 

B.  Reflected  Field  from  the  Strip 

The  reflected  fields  are  grouped  in  the  same  way  as  the  source 
fields  and  can  be  expressed  by 


FR(0) 


(4-3) 


where 


fS(0)  = 


M 

l 

m=l 


km 


•  Ay 


0 


for  ej,  <  e  <  e«2 
elsewhere. 


and 

(4-4) 


and 


*km  *  k^xkmcos0+ykms1n0^  ’ 
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r 


with  (x,|.n,yJin)  being  the  coordinates  of  the  image  point  for  the 

R  R 

source  point  (0^^)*  Note  that  0^  and  are  the  shadow  bour 
angles  of  the  reflected  field  for  the  kth  subaperture. 


Diffracted  Field 


As  previously  discussed,  the  incident  field  Ekj-  on  the  edge  QCj 
of  the  strip  scatterer  from  the  kth  subaperture  is  the  sum  of  all 
the  contributions  from  M  elements  in  this  subaperture.  That  is 


-^kmj 

A^km^  r- — - 

J  skmj 


(4-5) 


where  smkj.  represents  the  distance  from  the  source  point  (0,ykm)  to 
the  diffracting  edge  Q  •.  Note  that  j=l  or  2  indicates  the  two  edges 

^  J 

making  up  the  strip.  The  diffraction  coefficient  Dkj.  for  the  kth 
subaperture  may  be  expressed  by 


Dkj  "  Dl(skj'*kj’*kj)  +  ^^kj’^j^kj)  (4"6) 

where  <j>kj-,  <j>£j,  and 
of  the  subaperture  and  the  appropriate  edge,  as  shown  in  Figure  11. 
So,  the  pattern  of  the  diffracted  field  Fc*(e)  is  given  by 

Fdle)  ■  hJj(b)  •  DkJ  (4-7) 

where  is  defined  in  Equation  (3-12). 


skj-  are  determined  relative  to  the  midpoint  (0,yk) 
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CHAPTER  V 


TWO-POINT  GTD  SOURCE  FIELD  METHOD 

In  this  chapter,  two  types  of  source  distribution  are  considered 
and  analyzed  by  two  different  GTD  approaches.  The  first  one  is  a 
uniform  source  distribution  on  a  source  strip,  which  can  be  approxi¬ 
mated  by  the  diffracted  fields  from  two  edges  of  this  strip  illumi¬ 
nated  by  a  uniform  plane  wave.  The  second  one  is  a  source  distri¬ 
bution  with  zero  edge  illumination,  e.g.,  a  cosine  distribution,  which 
is  analyzed  in  terms  of  the  slope  diffractions  from  two  edges  of  the 
strip  with  the  same  illumination.  Both  approaches  are  more  efficient 
than  the  aperture  integration  method  as  will  be  shown  later. 

A.  Uniform  Source  Distribution 

To  start  with,  consider  an  antenna  source  with  a  uniform  distri¬ 
bution.  The  procedure  for  analysis  is  about  the  same  as  that  in  Chap¬ 
ter  III,  but  the  source  field  is  treated  by  a  two-point  GTD  approach 
as  described  below. 

1 .  Source  field  without  scatterer 
a.  Near  field 

The  near  field  of  an  aperture  antenna  as  shown  in  Figure  5  can  he 
analyzed  by  using  GTD.  In  which  case,  the  source  field  is  approximated 
by  the  diffraction  fields  from  the  two  edges  of  the  antenna  as  shown  in 
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figure  12.  Note  that  the  diffracted  fields  are  those  of  a  conducting 
strip  illuminated  by  an  incident  plane  wave. 

The  diffracted  field  from  each  edge  is  given  by 

-jksi 

E?(R,e)  =  Ei(Q  ,)  •  0n  •  e__ 

(5-1) 

"Jks2 

e|(r, e)  =  e1'(qs2)  •  oI2  • 

jj- 

where  Ei(Q$1)  =  Ei(Q$2)  =  ^ 

^11  =  ^1  (S1 » 4>i+<l)i )  (5-2) 

^*12  =  ( s2’ ^2+<^2^ 

which  are  defined  in  Equation  (2-16),  =  't'-j  =  and  v2,  Si 

and  s2  are  shown  in  Figure  12.  Note  that  only  d|  terms  are  used  here 
which  correspond  to  the  reflected  shadow  boundaries  in  the  z>0  half 
space.  It  has  been  shown  that  the  d|  terms  correspond  to  the 
physical  optics  or  aperture  integration  fields  [7], 
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0Si 

/  INCIDENT  PLANE  WAVE  E 

atL’ 

/ 

E'  (QS)) 

Figure  12. 

The  near  field  of  an  antenna  by  using 
two-point  GTD  approach. 

b.  Far  field 

In  the  case  of  R-*30,  as  shown  in  Figure  13,  the  far  field  approxi 

mat  ion  is  given  by 

E 


So  one  obtains 


Es(R,0)  =  Ei(Qsl)OI(s1;1r-0) 


-jks. 


-jks. 


-  +  E1  (Qso)Dt (Spiir+e) 

J*i  Jh 


(5-5) 


where  e’(Qs1)  =  E^Qc?)  =  En  in  this  case. 


%2'  o 


Then  the  far  field  pattern  function  Fs(0)  is  given  by 


F  (8)  =Eo 


jky  -I sin0  jky  ?sine"| 

DjfSjjir-eJe  +  DI(s2;Tr+6)e  J 


(5-6) 


as  Sp  s^-*®. 

2.  With  the  strip  scatterer  present 

The  geometry  of  the  antenna  in  the  presence  of  the  strip  scat¬ 
terer  is  shown  in  Figure  14.  Let  s^-  and  O^j  represent  the  distance 
and  the  angle  from  Q ^  to  the  scatterer  edge  Q  •  respectively,  where 
i=l  or  2  and  j=l  or  2.  It  follows  directly  that 

s1j  a7<xcj-St)2t^cj-ysi)2  ’  and  (5-7) 


(5-8) 
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Figure  13.  The  far  field  of  an  antenna  using 
two-point  GTD  approach. 


Figure  14.  The  geometry  of  the  two-point  GTD 
with  the  scatterer  present. 


a.  Source  field  with  scatterer  present 


Let  F^(e)  and  F^(0  )  be  the  source  field  patterns  from  each  an¬ 
tenna  edge  Qs^  and  in  the  presence  of  the  strip  scatterer,  re¬ 
spectively.  Similar  to  Equation  (3-4),  it  is  easy  to  see  from  Figure 
15  that 


F-(0) 


*0 


’°Iie 


jkysisin  0 


for  <  0  <  0^ 
elsewhere. 


(5-9) 


where  and  are  defined  in  Equation  (5-2).  By  superposition, 
the  total  source  field  pattern  Fs(0)  is  given  by 

Fs(0)  =  F*  (9)  +  f|(0).  (5-10) 

As  shown  in  Figure  16,  in  regions  I  and  V,  there  are  source  fields 
from  both  and  Q^.  In  regions  II  and  IV,  there  is  only  one  source 
field  from  either  Q$2  or  Qs^.  There  is  no  source  field  in  region 
III.  Note  that  there  may  be  no  region  III  in  the  far  field  if  the 
scatterer  is  small  compared  to  the  source. 

b.  Reflection  from  the  scatterer 

Then,  similar  to  Equation  (3-6),  the  contribution  of  the  reflected 
field  from  point  is  given  by 


,  and 


F-(0)  =i 


C  J'lfej 

+Eo°ii(R*1t'0R)  e 


for  9R2<  9  <  0^ 


elsewhere. 


(5-11) 


where  i  =  1  or  2, 

R  *  °°  for  the  far  field, 

0R  =  2AQ-0-1T , 

A0  is  the  tilt  angle  of  the  strip  scatterer 
and  <)>Ri  =  k(xIicos0+yIisin9). 

R  R 

Note  that  0^  and  are  the  shadow  angles  of  edge  Q^.  as  shown  in 
Figure  17.  The  reflected  field  regions  are  shown  in  Figure  18.  In 
regions  II  and  IV,  there  are  reflected  fields  from  the  edge  only. 

In  region  III,  there  are  reflected  fields  from  both  edges.  There 
are  no  reflected  fields  in  regions  I  and  V.  Again,  there  may  be  no 
region  III  in  the  far  zone. 

c.  Diffraction  from  the  scatterer 
Since  each  source  point  on  the  antenna  illuminates  both  edges 
of  the  scatterer,  there  are  4  diffracted  terms.  The  field  from  source 
point  Qsj  incident  on  the  edge  Qcj  of  the  scatterer  is  given  by 


Eo°si 


for  i=l,2  and  j=l,2. 


(5-12) 
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ANTENNA 


17.  The  reflected  field  of  the  two-point  GTD  method. 

IMAGE  OF  ANTENNA 
IN  STRIP 
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/ 

18.  The  illuminated  regions  of  the  reflected  fields. 


w 


1 


where  Oji  and  S^.  have  been  defined  in  Equation  (5-2)  and  Equation 
(5-7)  respectively.  The  pattern  function  of  the  total  diffracted 
field  Fd(0)  is  the  sum  of  the  4  terms  such  that 


2  2 


F°(e)  -  x  x  n>), 

1=1  j=i  j 


(5-13) 


where  F^(0)  corresponds  to  the  diffracted  field  of  source  from 


edge  Q  which  is  given  by 

L  J 


Fdj(0)  =  e\.  •  Dij  e^  for  1=1,2  and  j=l,2. 


(5-14) 


Similar  to  Equation  (3-10),  the  diffraction  coefficients  are  given  by 


D- .  =  Dt(L- -  D.(L .  .;<b.  .+6'. .) 
ij  Iv  ij’yij  yij'  r  ij’  ij  yij' 


(5-15) 


and 


Lij  = 


S1J*Sj 

stj*sj 


s . .  as  s 

ij  d  j 


(5-16) 


Note  that  <t>^ j  and  4>. ^  are  the  incident  and  diffracted  angles  from 
Qsj  to  Qcj,  respectively,  as  shown  in  Figure  19. 
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si  opr  diffraction  as  long  as  it  is  reasonably  far  from  the  main  beam. 
An  incident  plane  wave  E1  with  a  cosine  distribution  as  shown  in  Fig¬ 
ure  20  is  used  as  an  example  in  this  section,  where  E1  =  cos(j^)  and 
the  width  of  the  aperture  antenna  is  given  by  L.  Then  the  slope  dif¬ 
fraction  from  one  antenna  edge  is  given  by  Equation  (2-26),  such  that 


3E^(Q  .)  'JkS-j 

Ei  "  3n  Dpl^si,(^i+^i^  J~T.  for 


or  2, 


where 


3Ei(Qsl)  3[cos(f-)] 


an 


T(-y) 


,4 


7T^ 

L’ 


3E\QS2)  3 [cos  (^-)] 


3  n 


ay 


y4 


=  -  3L 
L’ 


(5-17) 


(5-18) 


^1  =  ^2  =  ?  ’  ^1=?”  + 


and 


3Dl(si  ;<t>  i 


(5-19) 


Note  that  s1  =  R  -  y^sine,  and  Sg  =  R-yS2sin9.  The  total  source 
field  Es  is  given  by 
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Figure  20.  A  cosine  source  distribution  modeled 
by  two-point  GTD. 


So  the  far  field  pattern  function  Fs(0  )  is  obtained  from  Equation 
(5-17)  as 


s  so  [  Jkyslsine 

F  (e  }  I  DpI(sl;iI_e)e  +  Dpi(s2;w+e)e 


jkys2sine 


(5-20) 


where  sQ 


9n 


which  is  equal  to  -  £  for  the  cosine  distribution. 


2.  Total  solution  with  the  strip  scatterer  present 
The  source,  reflected,  and  diffracted  fields  are  evaluated  in 
the  same  way  as  in  the  previous  procedure  except  that  the  source  fields 
correspond  to  those  of  the  slope  diffractions  from  the  antenna  edges. 
The  equations  for  the  various  field  contributions  are  given  in  the 
fol lowing  sections. 
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a.  Source  field  with  the  scatterer 


Let  F|(e)  represent  the  source  field  pattern  from  the  antenna 
edge  at  e^,  where  i=l  or  2,  as  shown  in  Figure  15.  Thus 


F-(0) 


0 

1 


for  9il<9<9(2 


,  and 


(5-21) 

s0  jkysisine 

3k  elsewhere. 


3  E 

where  sQ  =  g-^-  and  DpI  is  defined  in  Equation  (2-26),  and  all  other 
parameters  are  the  same  as  in  Equation  (5-9).  Then,  the  total  source 
field  with  the  scatterer  present  is  given  by 

F  (o)  =  l  f!(0).  (5-22) 

i=l  1 


b.  The  reflected  field 

Referring  to  Figure  17,  the  pattern  of  the  reflected  field  from 
the  fth  edge  of  the  source  is  given  by 

r 


F^e) 


x  S0 


%i 


+  DpI  ( R ;ir -e R )  e  for  ei2  <  e  <  e 


11' 


and 


(5-23) 


1.0 


elsewhere. 
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where  R.  and  tiR  are  the  same  as  those  in  Equation  (5-11).  The 
total  refle 
previously. 


total  reflected  field  F^(0)  is  the  sum  of  F^(0)  and  F^Q)  as  done 


c.  The  diffracted  field 


-I 


As  shown  in  Figure  19,  the  incident  field  £•■  is  given  by  Equation 

1  J 

(5-17),  such  that 


-I  so 


"jkSij 


Eij  =  31c  DPi(W+i> 


7s- 


(5-24) 


ij 


Then,  following  the  previous  procedure,  one  finds 


Ff j<  e>  *  e!j  -D>j 


(5-25) 


where  D-j  and  <|>dj  have  been  previously  defined  by  Equations  (5-15) 
and  (3-12).  Then  the  total  diffracted  field  pattern  is  given  by 


FSd(e>  ■  i  jj 


(5-26) 


CHAPTER  VI 


SLOPE  DIFFRACTION  FROM  THE  SCATTERER 

As  mentioned  in  Chapter  I,  when  the  strip  scatterer  is  located 
near  the  projected  aperture  of  the  antenna,  slope  diffraction  from 
the  edges  of  the  scatterer  must  also  be  included.  Large  aperture 
antennas  generally  have  rapid  field  variations  which  implies  that 
the  field  incident  on  the  edge  of  scatterer  may  likely  have  a  very 
significant  slope.  This  is  illustrated  in  Figure  34  in  the  next 
chapter.  It  was  found  that  slope  diffraction  models  the  rapid  source 
field  variations  very  well  for  a  uniform  source  distribution  if  the 
scatterer  is  not  too  close  to  the  projected  aperture.  However, 
slope  diffraction  does  not  work  so  well  for  a  cosine  source  distri¬ 
bution  (see  Figure  38)  because  the  source  field  varies  too  rapidly. 

The  geometry  used  to  evaluate  the  slope  of  the  incident  field 
is  shown  in  Figure  21,  where  the  latter  is  numerically  approximated 
by 
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with  Ei -,  Ei An,  0* and  si.  defined  in  Figure  21.  So  the  con- 
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tributions  from  the  slope  diffraction  terms  are  given  by 


T 


SCATTERED 


Figure  21.  The  geometry  used  to  calculate  the  slope 
of  the  incident  field. 
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where  the  slope  diffraction  coefficient  P-j  =  ^D.j.j/^<|>,  is  given  by 


pij  °pl ^Si ^i '  Dpl ^ si j ; ^i j+<t>i j ^  ’ 


(6-3) 


where  s-,  <(> i  and  ^ ^  are  shown  in  Figure  19,  and  DpJ  is  defined 
by  Equation  (2-26).  Recall  that  the  "+"  sign  is  used  for  the  TE  case 
and  the  sign  for  the  TM  case. 
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CHAPTER  VII 


RESULTS  AND  DISCUSSIONS 


'Ip  the  previous  chapters,  two  approaches  are  developed  for  cal¬ 
culating  the  effect  of  a  strip  scatterer  on  the  antenna  pattern: 
a)  the  large  subaperture  technique  and  b)  the  two-point  GTD  source 
field  method.  Both  approaches  use  the  GTD  method  to  calculate  the 
scattering  from  structures  in  a  manner  similar  to  that  used  in  the 
present  codes [  1,2].  In  this  chapter,  results  are  given  to  show  the 
effects  and  limitations  of  the  two  approaches.  The  results  of  the 
large  subaperture  technique  for  the  source  field  are  compared  in  Sec¬ 
tion  A  with  those  from  the  original  GTD  method  which  uses  super¬ 
position  from  small  elements.  The  accuracy  and  efficiency  of  the 
two-point  GTD  source  field  method  are  shown  in  Section  B.  Next,  when 
the  strip  scatterer  is  located  near  the  projected  aperture,  slope 
diffraction  can  be  used  to  improve  the  results  for  the  GTD  source 
field  methods  as  discussed  in  Chapter  VI.  The  effect  of  slope  dif¬ 
fraction  from  the  scatterer  and  its  limitatations  are  shown  in  Section 
C.  Finally,  the  effect  of  the  various  scatterer  geometries  and  lo¬ 
cations  are  shown  in  Section  D.  A  10X  source  antenna  is  used  in  these 
results  and  two  types  of  source  distributions  are  included:  type 
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1  is  for  the  cases  of  uniform  source  distribution  and  type  2  is  for 
the  cosine  source  distribution. 

A.  Large  Subaperture  Technique 

As  mentioned  in  Chapter  I,  the  large  subaperture  technique  works 
well  as  long  as  the  size  of  the  subaperture  is  no  larger  than  1  to 

2  wavelengths. 

1.  Type  1  -  Uniform  source  distribution 

The  source  field  of  the  uniform  distribution  (type  1)  without 
the  scatterer  present  is  shown  in  Figure  22.  Note  that  it  exists 
in  the  region  of  -  90°  <.  0  £90°  only.  The  far  field  patterns  are 
shown  in  Figures  23  and  24  for  two  different  scatterer  geometries. 
Both  cases  include  lA  and  2A  subaperture  sizes  and  are  compared  to 
the  results  using  small  A/3  elements.  Note  that  there  is  almost  no 
difference  between  the  results  using  A  subapertures  and  A/3  elements 
as  shown  in  Figures  23a  and  24a.  This  method  works  very  well  for 
most  regions  even  with  2A  subapertures  and  used  only  1/6  of  the  orig¬ 
inal  CPU  time  as  shown  in  Figures  23b  and  24b. 

2.  Type  2  -  Cosine  distribution 

A  cosine  distribution  is  used  to  represent  a  source  distribution 
with  zero  edge  illumination  is  these  examples.  The  results  for  type 
2  are  shown  in  Figures  25,  26  and  27  for  the  same  antenna/scatterer 
geometries  as  were  shown  for  type  1.  The  large  subaperture  technique 
does  not  work  so  well  as  it  does  for  type  1  because  the  source  field 
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pattern  of  a  10X  antenna  with  uniform  source 
(type  1  antenna). 


SCATTERED 


pattern  of  a  lOx  type  2  antenna  with  strip  scatter  B 


varies  more  sharply.  For  the  result  shown  in  Figure  26,  the  pattern 
value  at  angle  -151°  was  analyzed.  Table  I  compares  the  individual 
contributions  from  each  2A  section  of  the  aperture  (by  using  A/3  sub- 
apertures  which  is  equivalent  to  the  original  method)to  the  contri¬ 
butions  when  the  large  subaperture  method  is  used  with  2A  subaper¬ 
tures.  Although  the  agreement  between  the  contributions  from  each 
section  is  good,  there  is  a  large  difference  in  the  total  fields. 

Note  that  the  total  field  is  small  compared  to  the  individual  con¬ 
tribution  as  seen  from  Table  I.  This  cancellation  effect  requires 
a  very  high  accuracy  for  the  individual  contribution  in  order  that 
the  total  field  be  accurately  calculated. 

B.  Two-point  GTD  Source  Field  Method 

The  source  field  patterns  without  the  scatterer  present,  as  ob¬ 
tained  by  using  two-point  GTD,  are  shown  in  Figures  28  and  31  for 
source  types  1  and  2,  respectively.  Note  that:  (1)  GTD  breaks  down 
near  the  main  beam  and  (2)  the  source  fields  in  the  back  ( i . e . ,  for 
|o|  >  90°)  are  omitted  for  the  2-point  GTD  in  order  to  compare  with 
the  aperture  integration  source. 

Figures  29,  30,  32,  and  33  show  the  results  of  using  the  two- 
point  GTD  method.  These  results  are  compared  to  the  results  of  using 
aperture  integration  (AI)  for  the  source  fields.  It  is  apparent  that 
2-point  GTD  for  the  source  field  works  very  well  in  almost  all  regions 
except  the  main  beam,  where  the  GTD  breaks  down.  The  efficiency  of 
the  2-point  GTD  method  can  be  obtained  by  comparing  the  computer  CPU 
time  as  follows.  By  using  the  AI  source  with  A/3  elements  for  a  10A 
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Figure  30.  The  far  field  pattern  of  a  10*  type  1  antenna  with  scatterer 


ie  far  field  pattern  of  a  lOx  type  2  antenna  with  scatterer 


diiluimd  (30  elements),  it  takes  about  110  seconds  to  get  the  result 
shown  in  Figure  29 <  Since  the  CPU  time  is  about  proportional  to  the 
number  of  elements  in  the  antenna,  it  will  take  more  time  if  the 
antenna  becomes  larger.  But  by  using  two-point  GTD  method,  no  matter 
how  large  the  antenna  is,  it  only  takes  about  12  seconds. 

C.  Slope  Diffraction  from  the  Scatterer 

As  discussed  in  Chapter  VI,  slope  diffraction  from  the  scatterer 
is  significant  in  some  regions  when  the  GTD  source  method  is  used 
and  the  scatterer  is  near  the  projected  aperture  of  antenna.  For 
type  1,  if  one  edge  of  the  scatterer  is  located  3A  above  the  projected 
aperture  of  the  antenna,  slope  diffraction  can  almost  perfectly  elimi¬ 
nate  the  slope  discontinuity  at  the  shadow  boundary  near  e=-110°  for 
the  reflected  field  as  shown  in  Figure  34.  But  when  the  scatterer 
is  located  only  IX  or  2X  above  the  projected  aperture,  slope  diffrac¬ 
tion  is  not  sufficient  to  solve  the  problem,  as  shown  in  Figures  35 
and  36.  The  reason  is  that  the  source  field  is  varying  too  sharply 
near  the  shadow  boundary  and  can  not  be  represented  by  a  uniform  plus 
slope  field.  For  a  cosine  distribution,  the  incident  field  on  the 
edges  of  the  scatterer  varies  even  more  sharply  and  consequently, 
the  problem  becomes  even  more  serious.  As  shown  in  Figures  37  and 
38,  slope  diffraction  is  not  enough  to  solve  the  problem  in  the  small 
region  near  the  reflected  shadow  boundaries  until  the  scatterer  is 
located  about  6A  above  the  projected  aperture  for  this  case. 
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(9,12) 


e  diffraction  from  the  scatterer  E  for  type  2  antenna 


Combination  of  Large  Subaperture  and 
Two-point  6TD  Source  Methods 


In  this  section,  the  results  of  the  previous  sections  are  pre¬ 
sented  to  show  the  effect  of  the  scatterer.  The  results  were  obtained 
by  using  two-point  GTO  method  in  all  regions  except  the  main  beam, 
where  the  large  subaperture  method  is  used.  In  order  to  see  the  ef¬ 
fect  of  the  scattering,  the  source  fields  without  the  scatterer  present 
are  plotted  as  dashed  lines  in  Figures  40  to  46.  Note  that  the  back 
source  fields  for  | e|  >  90°  of  the  2-point  GTD  are  kept  in  this  sec¬ 
tion  as  shown  in  Figures  39  and  42.  It  is  clear  that  each  scatterer 
has  its  maximum  effects  near  the  reflected  and  blocked  regions,  which 
are  marked  on  the  figures. 
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LO.  type  1  antenna  with  scatterer  A.  Two  approaches 
compared  to  the  source  field. 


CHAPTER  VIII 


EXAMPLES  OF  OTHER  APPLICATIONS 

In  Chapters  V  and  VI,  only  a  normally  incident  plane  wave  is 
considered;  however,  this  is  not  a  necessary  condition.  In  this  chap¬ 
ter,  some  more  general  cases  are  considered  and  two  examples  are 
demonstrated  as  follows: 

A.  Scanning  Mode 

The  first  example  is  a  plane  wave  incident  at  an  angle  0m  for 
a  phase  scanned  planar  array  as  shown  in  Figure  47.  The  diffraction 
solution  is  modified  in  terms  of  the  following  angles: 

*i  a  7  "  em  »  4>2  =  7  +  em’^l  =  7  -e»  and  ^2  =  7  +  0-  (8_1) 

Considering  only  the  DI+  terms  as  done  previously,  one  finds 

4>1  +  4»i  =  it  '  (0  +  9m),  and 

(8-2) 


*2  +  ^  "  +  (a  +  em), 

while  the  remaining  terms  are  the  same  as  those  shown  in  Figure  48. 
From  Equation  (5-6),  as  s^,  s^  -*■  the  pattern  function  Fs(e)  of 
the  source  field  is  given  by 
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c  jky0sine 

Fs(e)  =  Eo[DI(s1,tr-{64em))e  SA 


jky  ?sin0 

+  DI(s2;Tr+(0+em))e  ] 


(8-3) 

Note  that  this  equation  is  obtained  from  Equation  (5-6)  by  changing 
0  to  e+om. 

With  the  presence  of  the  strip  scatterer  as  shown  in  Figure  48, 
the  source,  reflected,  and  diffracted  fields  can  be  calculated  as  in 
section  B  of  Chapter  V  except  that  the  angles  associated  with  the 
diffraction  coefficients  should  be  changed  from  8  to  0+em,  where 
0m  is  the  scan  angle. 

B.  Cylindrical  Wave  Radiator 

The  second  example  is  a  line  source  behind  a  conducting  slot 
as  shown  in  Figure  49a.  The  analysis  is  similar  to  that  used  to 
study  horn  antennas  as  done  in  [8].  This  problem  can  be  mod¬ 
eled  by  a  line  source  located  at  (xs,ys)  in  the  presence  of  a  per¬ 
fectly  conducting  strip  as  shown  in  Figure  49b.  In  this  case,  the 
source  fields  of  the  strip  scatterer  include  two  contributions,  one 
is  the  geometric  optics  term  Eg  q  and  the  other  is  the  diffractions 
from  edges  Qs^  and  Qs2  as  discussed  previously.  Then  the  total  field 
is  the  superposition  of  the  fields  caused  by  these  two  contributions 
as  discussed  below. 
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1.  Terms  Caused  by  Geometric  Optics  Source 

First,  the  geometric  optics  field  without  the  strip  scatterer 


present  is  given  by 
ejks' 

EG*°*  -7^ 


V  o 


for  8m  <  0  <  0m2,  and 
elsewhere. 


(8-4) 


where  Om2,0m2>  and  s‘  are  shown  in  Figure  50.  Then,  with  the  strip 

scatterer  present,  the  source,  reflected,  and  diffracted  fields  caused 

r  o 

by  E  can  be  calculated  separately  as  before.  Note  that  depending 
on  the  source  point  (xs,y$)  and  the  location  of  the  scatterer,  some 
of  these  contributions  may  not  exist.  A  case  in  which  all  the  con¬ 
tributions  exist  is  shown  in  Figure  51a.  On  the  other  hand,  only 
the  source  field  and  part  of  the  reflected  and  diffracted  exist  for 
the  geometry  shown  in  Figure  51b.  Finally,  an  example  without  any 
reflected  and  diffracted  fields  caused  by  the  geometric  optics  source 
is  shown  in  Figure  51c. 

2.  Terms  Caused  by  Diffraction 
Sources 

Next,  the  contributions  due  to  the  source  fields  caused  by  the 
diffractions  from  and  QS2  can  be  analyzed  as  before  except  that 
<j>j,  <t>2«  E^Q^)  and  E^Q^)  are  appropriately  modified.  The  geometry 
associated  with  this  case  is  shown  in  Figure  52,  where  4>j,  $£,  sj 
and  s£  are  determined  by 
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Besides,  as  defined  before,  ^  -  01  and  <j>2  =  +  @2.  So  the  source 

field  Ep(r,0)  at  point  P  due  to  the  diffraction  from  Q$1-  given  by 


-jks! 


e_ 


-jksi 


E  (r,0)  =  — —  DI(Li;«|)i+4,i ) 

/c  /si 


for  i  =  1,2,  and 


si 


(8-7) 
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If  one  considers  the  far  field  of  the  source  field,  that  is, 
when  r-K*.,  then  s^,  s2-*»  and  0^=02=0.  The  pattern  of  the  source  field 
becomes 


Fs(e) 


-jks{ 

of _ 


jky  ,sino 
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jky  2sine 

^  I  ^  s2  ’  <^2+<^>s ) e 


(8-8) 


Note  that  it  is  similar  to  Equation  (5-6)  except  that  eVQ^)  is  changed 


CHAPTER  IX 


CONCLUSIONS 


As  discussed  in  Chapter  III,  GTD  can  be  used  to  accurately  an¬ 
alyze  scattering  problems.  However,  the  use  of  the  aperture  inte¬ 
gration  source  is  not  efficient,  especially  when  the  antenna  is  large 
in  terms  of  the  wavelength.  The  large  subaperture  technique  and  two- 
point  GTD  source  field  method  were  studied  in  this  thesis  as  ways 
to  improve  the  efficiency,  and  their  corresponding  advantages  and 
limitations  are  shown  in  the  results.  From  these  results,  it  can 
be  concluded  that  the  subaperture  technique  works  well  for  lx  sub¬ 
apertures  and  that  it  works  reasonably  well  for  2A  subapertures  except 
in  some  small  regions  for  certain  cases.  Two-point  GTD  turned  out 
to  be  a  very  good  approach  for  both  types  of  aperture  distributions 
(uniform  and  cosine)  in  all  regions  except  the  main  beam  as  long  as 
the  scatterer  is  not  too  close  to  the  projected  aperture  of  the  an¬ 
tenna 

If  the  scatterer  is  near  the  projected  aperture,  the  slope  of 
the  incident  field  on  the  edges  of  the  scatterer  becomes  significant. 
In  this  case  slope  diffractions  from  the  scatterer  can  be  added  to 
improve  the  accuracy  near  '■  shadow  boundaries,  but  it  can  not  solve 
the  problem  completely  it  i..ie  scatterer  is  too  close  to  the  projected 
aperture  of  the  antenna. 
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Since  the  purpose  for  developing  the  two  approaches  discussed 
in  this  thesis  was  to  improve  the  efficiency  of  the  calculations, 
consequently  the  CPU  time  is  another  important  result  to  discuss. 

For  the  large  subaperture  technique,  the  efficiency  improvement  is 
approximately  proportional  to  the  number  of  the  elements  used  in  one 
subaperture.  That  is,  if  3  elements  are  used  in  one  subaperture, 
it  only  takes  about  1/3  of  the  original  CPU  time  to  calculate  the 
results.  For  the  two-point  GTO  method,  the  CPU  time  is  independent 
of  the  size  of  the  antenna.  For  the  results  shown  from  Figures  29 
to  33,  the  GTD  method  is  more  efficient  than  using  A/3  elements  by 
a  factor  of  about  9  for  a  10A  source  antenna.  It  is  obvious  that 
the  efficiency  will  be  even  better  for  larger  source  antennas. 

In  this  thesis,  only  a  single  2-D  strip  scatterer  was  studied. 
If  3-D  cases  or  multiple  strip  scatterers  are  considered,  more 
complicated  factors  such  as  corner  diffraction,  polarization,  second 
order  diffraction,  etc.  should  also  be  considered. 
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APPENDIX  A 


ELEMENT  PATTERN  USED  FOR  APERTURE  INTEGRATION 


As  discussed  in  Chapter  V,  the  radiation  pattern  of  an  aperture 
antenna  can  be  analyzed  by  using  the  2-point  GTD  method.  In  this 
section,  a  modified  Huygen's  element  pattern  cos  e/2  is  derived  by 
comparing  the  results  from  2-point  GTD  method  and  from  aperture  in¬ 
tegration  (AI).  Using  this  element  pattern  in  the  AI,  one  gets  a  good 
agreement  between  GTD  and  AI.  By  using  GTD  as  shown  in  Figure  A-l, 
the  radiation  pattern  EqTD  is  given  by 


EGTD 


Dju-e) 


-jkSj  -jk$2  ■ 

*  0j(,+8) 

Jh  Jh  J 


=  E 
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where 


D|(n±rt) 


F(kLa) 
sin  0/2 


and  k  is  wave  number  =  2w/X 


(A-l) 


92 


Since  F(kLa)  -►  1  for  the  far  field, 


E1  e"J  * 


4n  sin  | 


—  [a  sin  e] 


_  Eie  j.  sin (ip  sin  e) 


sin  | 


(A-2) 


On  the  other  hand,  by  using  aperture  integration  as  shown  in  Figure 
A-2  and  representing  the  element  pattern  as  F^( e) ,  we  obtain 


?  -jks' 

eai  =  c  4  J.  ~f=~  w 


where  C  is  a  constant.  For  the  far  field,  s'  =  R-y'-sin  e 


eJky  sine  2  g-jkR 


(A-3) 
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(A-4) 


Using  the  Huygen's  source  element  pattern  as  given  by 
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Fe<6)  -  ^ 


results  in 


E 


AI 


e-jkR 


However,  simple  use  of  a  modified  Huygen's  source 
FE(e)  =  cos  | 


gives  a  result  which  is  identical  to  GTD  as  follows: 


(A-5) 


(A-6) 


(A-8) 
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